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1 Introduction 

t*"- | Practical calculation of transport properties of quantum networks is often 

reduced to the scattering problem for a one-dimensional differential operator 
^ ■ on a quantum graph, see for instance [H El [S] . Quantum graph plays a 

role of a solvable model for a two-dimensional network, see [HI El IB]- Basic 
detail of the model is a star-shape element with a self-adjoint boundary 
condition at the node. It was commonly expected that the realistic boundary 
condition is defined by the angles between the wires at the node. For instance 
| the boundary conditions for the T-junction, pQ, is presented in terms of limit 

values of the wave-function on the wires {^if-i := ip and the values of the 
■ corresponding outward derivative {^}? =1 := ip' at the node: 

O - 

^t- P~ Vi = Va = ^ , M + V4 + V4 = o, (i) 

o . 

or in the form 
Q-T with the projection 



> 

The scattering matrix of such a junction is S = I — 2Pq, see [H HI [TO]. In 
PQ (5 is interpreted as a free parameter describing the strength of coupling 
between the leg and the bar of the T-junction. In [11] the condition ( [2]) is 
used for analysis of spin-dependent transmission across the quantum ring. 

We derive the boundary condition ([2]) and similar condition for any junc- 
tion from [13] and interpret the corresponding free parameters P,Pq. 



P V=0, P ^ = (2) 



p 2 p 

p 1 1 
p 1 1 



2 Intermediate Hamiltonian 

Consider one body scattering problem on the junction £1 formed by few 2-d 
equivalent semi-infinite wires fL-, s = 1, 2, . . ., attached to the quantum well 
£lint via the orthogonal bottom sections T s . The corresponding one-body 
Hamiltonian for the spin- less electron is scaled, via replacement of energy 
E by the spectral parameter A = h~ 2 2m * E, to the standard Schrodinger 
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operator with zero conditions on the boundary d£l, and a constant potential 
in the wires q\n s = qoo- 

Ci/i = - A ip + qtp. (3) 

We assume, following [13], that the potential on the well is defined by the 
scaled constant electric field £: q\n int = Qint(x) = (£,§),§ S tg)\\ u . The role 
of the non-perturbed Hamiltonian is played by the Schrodinger operator 
C ou t on the wires with zero boundary conditions on the union V = L) S T S of 
the bottom sections, which play roles of solid walls, separating the well from 
the wires: 

V'lr = 0. (4) 

The eigenfunctions of £ ou t in the wires S7 S : < x s = x < oo, < y s = y < 5 
are combined of running waves 

exp(±i\/ A — ir 2 l 2 5~ 2 — qoox) y/2/5 sin irly/8 

:= exp(±ipix) e s . 

Hereafter we use on the wires the local coordinates x,y, < x, oo, < y < 
8. The eigenfunctions of the operator Li nt defined by (|3|4p on the quantum 
well Flint are standing waves. Replacement of the solid wall condition @ by 
the matching condition is a strong perturbation blending the standing waves 
on the quantum well with the running waves in the wires. This is a pertur- 
bation on the continuous spectrum, so the convergence of the corresponding 
series can't be estimated in spectral terms of self-adjoint operators. In [13] 
we suggested a modified analytic perturbation procedure based on introduc- 
tion of an Intermediate Hamiltonian obtained via appropriate splitting, see 

[E] of c. 

Assume that the Fermi level in the vires lies in the middle of the first 
spectral band 2m* Eph~ 2 := \f = ir 2 5~ 2 5/2 + q 00 . In that case all branches 
of the continuous spectrum with thresholds n 2 l 2 5~ 2 + q^, I > 2, are closed, 
that is all exponential solutions exp±ip/x e^, , I > 2, of the homogeneous 
equation —Aif^+qip = Xip in the wires are exponentially decreasing. Impose, 
additionally to (pQ) the semi-transparent boundary condition on the bottom 
section 

3 

u\ r ±E + , E + = \/ e], (5) 

s=l 

This boundary condition prevents excitations from the first channel 7i+ = 
E + x Li2{R+) in the wires from entering into the quantum well, and, vice 
versa, exiting from the quantum well into the wires, but does not stop the 
excitations from the closed channels 71- = [Vz>i e i] x L2{R+). The corre- 
sponding operator C F is split into the orthogonal sum of the trivial operator 
in the open first channel 

U = ~dx 2 + + n ' 

with zero boundary condition on T, and the intermediate Hamiltonian L F 
in the orthogonal complement TL F = {Z/2(f2) Tt+}. The spectrum a a (l F ) 
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Figure 1: Spectrum of L F 

coincides with the first spectral branch and the continuous spectrum a a (l F ) 
of L F is [47r 2 /<5 2 + Qoo, oo). There is a finite number of eigenvalues of L 
situated inside the first spectral band Ai = [tt 2 /S 2 + (/oo, 47r 2 /<5 2 + qoo] and 
a countable number of embedded eigenvalues of X F £ o a (L F ) accumulating 
at infinity. 

3 Scattering matrix via Intermediate DN map 

Usually the Scattering matrix on the quantum network is obtained via 
matching exponential solutions in the wires to the solutions of the homoge- 
neous Schrodinger equation inside the vertex domain, see, for instance |16j . 
This approach requires solving an infinite algebraic system. We consider the 
boundary problem for the intermediate equation: 

- A i> + qip = Xip, ip G L 2 (lo), 

rP\r = Vr € E+, 3A ^ 0. (6) 

The solution ijj exists for all complex A and has normal limit values on the 
continuous spectrum. We introduce the Dirichlet-to-Neumann map of the 
intermediate Hamiltonian (DN-map) as 

VM F :#^P+^Vr, (7) 

with P + := Pe + - It is a 3 x 3 matrix- function which is obtained via differen- 
tiation with respect to exterior normal of the resolvent of the intermediate 
operator restricted onto £li n t and framed by P+. It has the kernel: 

It has the spectral representation on the complement of the spectrum of 
a(L F ) 

P dMZ) \ ip 

VN f {\) = Y j + 9n /V 9n +JC F (\), (8) 
— ' A — \ r 

r 

where the summation is extended over discrete spectrum of L F and K, F {\) 
contains an integral over the continuous spectrum of L F . The scattering 
matrix of C is obtained via matching of the scattering Ansatz on the open 
channel in wires with p = p\ : 

e Pi e + +e * P S(p)e + (9) 
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to the limit values on the spectrum, 9A — > of the solution of the above 
intermediate boundary problem (|7|): 

ip [e + - S(p)e + ] = VM F {\) [e + + S(p)e + ] . 



Solving this equation we obtain, see |13j . the formula for the scattering ma- 
trix of the operator C on the first spectral band A = [n 2 /5 2 +qi n fty, 4-7T 2 /5 2 + 
Qinfty, ] in terms of C F by the formula 

_ M'w + ip p + 
yl -VM F (\) + ipP + 

The DN-map T>J\f F of the intermediate Hamiltonian L F is connected with 
the standard [12J DN-map A of the operator Li nt on the quantum well by 
the formula 

VM F = P+VNP+ - P + VMP_ T>- l P^VNP+. (11) 

Here 

D = P^VMP- - K_ , 

with 

K- = — ®J2piPi, \K-\ > vr^ 1 vW^- 
l>l 

and P- := © Sz>i ^2 S e s)( e <r Near the eigenvalue Ao of Li nt the DN-map 
can be decomposed as 

VM F {\) = M% + y M% := M% + ^ (12) 

s^O 

where S = P_^S^-\r, s = 0,1,... are the boundary currents of eigenfunc- 
tions $ s , s = 0,1,2... of the operator L« n t, A s - the corresponding eigen- 
values. Spacing between the eigenvalues of Li nt is connected to the di- 
ameter di n t of Qint as po = min^o | Ao — K\ = 0(d~ 2 t . Due to the spec- 
tral estimate | A^o (Ao) | ~ Po 1 - Fot relatively thin networks the analytic 
perturbation procedure can be developed based on (fTT|) . with the small 
parameter S/di n t, to obtain VM F . Denote by ]C± t ± the matrix elements 
P±K,qP± with respect to the orthogonal decomposition £2 CO = E + © £L. 
If only one eigenvalue Ao of Li nt is situated of the essential spectral interval 
At = [X F — 2h 2 m*KT, X F + 2h 2 m*KT, we can represent the denominator D 
of (112p near Ao, with P-(f>o := 0g with a controllable error k: 

D(A) = + - K_ := ^% + k 

A — Aq A — A 



so that the whole expression (jlip can be calculated via analytic perturbation 
procedure, since k « 5 _1 [1 + >> 1: 



p </»f ) (^0 



A (13) 



v 
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with Xq = Xq + (P_0o, k P-(j>o) ~ Ao and cp F = P+(/>o — Ah A; -1 P_<^o ~ 

P~\-(j)Q. For low temperature only electrons with energy close to Fermi level 
Ep contribute to transport phenomena. Hence T>Af F may be substituted 
by the single resonance term jMaJjjMa ; thus resulting in the approximate 
expression for the scattering matrix on At : 



ip p _ p o<t>p) Wo 
i P P + Ai i — ,\ 



Sapprox 

(14) 



'"F- 1 U ~i 

see [13 1 and more details in [17J. 



4 Boundary condition at the junction 

The approximate scattering matrix can be obtained from the energy-dependent 
boundary condition at the vertex imposed onto the scattering Ansatz (J9]) in 
the wires: 

ip[I - S approx {X)]ii) = [I + S a p prox (X)}ip' . (15) 

The polar terms in the numerator and in the denominator of (|14|) have the 
dimension cm" 1 and can be represented via the relevant one-dimensional 
orthogonal projection P := e ) (e with e := (ej, eg, . . . eg) =|| O 
4> := a _1 (/> . Then <p Q ) (<fi = a P . Denoting by P Q the complementary 
projection / — P Q in £2 CO, we obtain 



^approx (A) — -P Q + 



ip(A - A ) + a 
ip(X — Aq) — a 



P (16) 



^0 +©(A)P - 



The factor is close to —1 on the essential spectral interval At, for low 
temperature 2nm * \/3/2kT « bo? f? . Then, in first approximation, the 
energy-dependent boundary condition (fTST) is reduced on A7 1 to iP^ip — 
Potp' as 0, or, due to orthogonality of Po,P{j~, to Prj~ip ~ 0; -P ^' ~ 0. This 
condition coincides with the above Datta condition (pQ) presented in form 
d2]). Our analysis reveals the meaning of the projection Pq: it coincides 
with the projection onto the one-dimensional subspace defined by the vector 
P + |r of boundary values of the normal derivatives of the resonance 
eigenf unction, projected onto E + . This formula is valid also for general 
junctions. 



5 Example 

Consider a two-dimensional quantum network constructed as a simplest 
asymmetric T-junction of three straight quantum wires width ir/2 attached 
to the quantum well - the square Oj n t on £-plane : < £ x < 7r, < £ 2 < tt. 
Assume that the first wire uj\ = {—00 < £ t < 0, tt/2 < £ 2 < ir} is attached 
orthogonally to the left side of <9f2 on Ti = {£1 = 0, tt/2 < £2 < vr}, the 
second wire uj 2 = {0<£i<7r/2,7r<£2<°o}is attached in the middle of 
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Figure 2: Simplest asymmetric T-junction 

the upper side, and the third wire u> 3 = {ir < ^ < oo, < £ 2 < it} is at- 
tached to the middle of the right side of Oj n j. On the constructed quantum 
network 0, = Vti n t U wi U W2 U W3 consider the scattering problem for Lapla- 
cian with homogeneous Dirichlet boundary condition at the boundary. The 
cross-section eigenfunctions in the first channel in the wires uil, uj — 2, uj^ 
are : 



— sin 2£ 2 , 

\/7T 



r 2 



-^cos 2^i, 

\/7T 



—= COS 2£ 2 



i'l 



r 3 

The Dirichlet Laplacian on Qi n t has on the first spectral band Ai = [4, 16] 
the eigenvalues A = 5 ,Ai = 8, A 2 = 10 and A 3 = 13 with eigenfunctions 
$0) < J ) 2- C &3i ^0 = ~ sm £i sin 2^2- The boundary currents of 3>o are 



dn 



ri 



2 . 

— sin 2^2, 
7r on 



r 2 



— sin(i, 

7T 



9n 



7T 



• sin 2^2- 



Assume that the Fermi level of the material is situated between the first and 
the second thresholds of the network 4 < E F < 16 close to the eigenvalue 
Ao = 5 cm , X F = 4.33 cm~ 2 . The electrons are supplied to the network in 
the first spectral band from the second wire across the bottom section T2 and 
exit across Ti, T3. Due to orthogonality of the cross-section eigenfunction of 
the open channel to the boundary currents of the eigenfunctions <J> , $3 the 
corresponding modes are not excited. An essential link to the closed channels 
is supplied only by $o> the contribution from other eigenfunctions either 
vanish or are suppressed due to the factors (Ao — X s ) i n the denominator. 
The link of <3?o only to the closed channel in oj^ gives a scalar equation for 
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the re-normalized eigenvalue Ao of the intermediate Hamiltonian, since only 
the contribution from T 3 is non-trivial: 



d® 
p u o 



dn 



dT 3 +V16-5 :=D_(A ) 



11 [0.67 + A - 5] 



(17) 



From this equation we obtain the resonance eigenvalue of the intermediate 
Hamiltonian: Ao = 4.33 = X F . The boundary current of the corresponding 
eigenfunction essentially coincides with the boundary current of the nor- 
malized resonance eigenfunction $ of the Dirichlet Laplacian on f^o- The 
projections of the resonance boundary currents onto E + are 



dn 
d$ 



r 2 



dn 



r 2 



#2 



--^= = -0.56, 

'7T 



0.43, 



r d$ 


3 


lr 3 dn 


r 3 



16^2 
3^r^/^r 



0. 



Then the normalized vector of the boundary current is eo = (—0.8, 0.6, 0), 
and the boundary conditions at the junction for low temperatures are rep- 
resented by the formulae ([2]) with Po = e o) (e , which is different from the 
condition for a symmetric junction suggested in [TJ for symmetric T-junction. 
For the higher temperatures the boundary condition is energy dependent and 
can be represented in form (|15p . with the approximate scattering matrix 



Sappr \P) 



iV\ - 4P + 



0.15 



Po 



A-4.33 



iy/X - 4P + + 0.15 



A-4.33 



with P = e )(e . 
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